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Abstract: We have produced high transmission sub-wavelength tapered 
optical fibers for the purpose of whispering gallery mode coupling in fused 
silica microcavities at 780 nm. A detailed analysis of the fiber transmittance 
evolution during tapering is demonstrated to reflect precisely the mode 
coupling and cutoff in the fiber. This allows to control the final size, the 
number of guided modes and their effective index. These results are checked 
by evanescent wave mapping measurements on the resulting taper. 

© 2008 Optical Society of America 

OCIS codes: (060.0060) Fiber optics; (060.2270) Fiber characterization; (140.3948) Micro- 
cavity devices; (140.4780) Optical resonators 



References and links 

1. F. Goiithier, S. Lacroix, X. Daxhelet, R. J. Black, and J. Bures, "Broad-band all-fiber filters for wavelength 
division multiplexing application," Appl. Phys. Lett. 54, 1290-1292 (1989). 

2. T. A. Birks, W. J. Wadsworth, and P. S. Russell, "Supercontinuum generation in tapered fibers," Opt. Lett. 25, 
1415-1417 (2000). 

3. J. C. Knight, G. Cheung, F. Jacques, and T. A. Birks, "Phase-matched excitation of whispering-gallery-mode 
resonances by a fiber taper," Opt. Lett. 22(15), 1129-1131 (1997). 

4. S. M. Spillane, T. J. Kippenberg, O. J. Painter, and K. J. Vahala, "Ideality in a fiber-taper-coupled microresonator 
system for application to cavity quantum electrodynamics," Phys. Rev. Lett. 91(4), 043,902 (2003). 

5. A. W. Snyder, "Coupling of modes on a tapered dielectric cyHnder," IEEE Trans. Microwave Theory Tech. 
MT18(7), 383 (1970). 

6. J. Love, W. Henry, W. Stewart, R. Black, S. Lacroix, and F. Gonthier, "Tapered Single-Mode Fibres and Devices 
Part 1 : Adiabaticity Criteria," lEE PROCEEDINGS 138, 343-354 (1991). 

7. P. N. Moar, S. T. Huntington, J. Katsifolis, L. W. Cahill, A. Roberts, and K. A. Nugent, "Fabrication, Modeling, 
and Direct Evanescent Field Measurement of Tapered Optical Fiber Sensors," J. Appl. Phys. 85, 3395-3398 
(1999). 

8. T. A. Birks and Y. W. Li, "The Shape of Fiber Tapers," J. Lightwave Technol. 10, 432^38 (1992). 

9. A. W. Snyder and J. D. Love, Optical Waveguide Theory (Chapman and Hall, London, 1983). 

10. "Spectrogram," Wikipedia URL h ttp : //en . wikipedia ■ org/wiki/Spectrogram| 

11. J. C. Knight, N. Dubreuil, V. Sandoghdar, J. Hare, V. Lefevre-Seguin, J. M. Raimond, and S. Haroche, "Charac- 
terizing whispering-gallery modes in microspheres by direct observation of the optical standing wave in the near 
field," Opt. Lett. 21, 698 (1996). 



Introduction 

Tapered optical fibers have attracted a large interest in the two last decades, due to their nu- 
merous applications as directional fiber- fiber couplers, filters for wavelength multiplexing |[T1, 
biosensors and more recently for the production of broad-spectrum-light from femtosecond 



laser impulsions IJ]- They have also been used successfully to selectively excite the whispering 
gallery modes (WGM) of spheroidal or toroidal microcavities The need for low losses 

has often led to work with waist diameters in the range 5 to 50 ^m while the weakness of the 
evanescent wave was, if needed, compensated by long interaction lengths. Efficient coupling to 
microcavities WGMs implies interaction lengths reduced to a few micrometers which in turn 
requires to work with smaller taper diameters and to control the effective index of the taper 
modes. The mode coupling and the resulting energy transfer have been extensively analyzed in 
the literature ||5]|6]|7l, either by Beam Propagation Method or, for slowly varying tapers radii, 
applying the coupled mode theory. The latter approach is best suited for tapers with character- 
istic lengths long enough to provide a physical meaning to the so called "local modes", defined 
at each position as the modes guided by a cylindrical fiber with the same diameter 

In this letter we describe how we produced sub-wavelength nearly adiabatic tapers with an 
overall transmittance higher than 90%, thanks to very efficient adiabatic transfer from the sin- 
gle mode of the untapered fiber to the fundamental mode of the central part of the taper In the 
frame of "local modes" modelling, we show that the final properties of the taper can be precisely 
deduced from the transmittance evolution during the tapering process. Moreover an experimen- 
tal criterion for a single-mode taper is demonstrated. Near-field mapping measurements of the 
tapers obtained by this method confirm our model. 

Taper fabrication 



Fig. 1. Top: Sketch of the experiment, showing the microtorch in the center, the two fiber- 
clamps FC, the rod R, the motorized translation stages MTS, the photodiode PD. Bottom: 
Shape of the resulting taper (not to scale), with definition of abscissa z, lengthening L, 
hot-zone h, taper waist w, initial radius tq. 

Our experimental setup is sketched in Fig.[T] The single mode FS-SN4224 fiber is held by 
two fiberclamps, sliding along two parallel stainless steel rods. The fiber is heated to its soft- 
ening point by a butane microtorch, using a specially designed nozzle producing a short flame 
which is « 10 mm wide along the fiber axis. The fiber position in the flame is a critical param- 
eter which can be adjusted by moving the nozzle using a three axis translation stage. The two 
fiberclamps are then symmetrically moved apart at the same velocity v of about 40 jtim • s^' by 
two motorized translation stages. During the whole process, the fiber transmission is monitored 
using a laser at the working wavelength X w 775 nm and a photodiode (PD). 

Figure [2] shows a typical recording of the fiber transmission during the tapering process, 
which lasts about 6 minutes, as a function of the fiber lengthening L = 2vt. After a slight de- 
crease of the transmission, oscillations appear, the frequency of which increases with time, as 
shown in the left inset. These oscillations are known to result from the interference between 
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Fig. 2. Fiber transmittance as a function of the fiber lengthening L (right scale). Note the 
95% final transmission (ie —0.22 dB insertion loss). The curve in solid squares is the ob- 
served decrease of the fiber waist vv (left scale). Left inset: zoom on the transmittance curve 
over 100 /xm for L = 20 mm and L = 30 mm. Right inset: zoom on the last amplitude drop. 



different modes and the modulation of their envelope is due to the beating of different frequen- 
cies. The most important feature is the amplitude drop occurring (for this taper) at L = 35 mm. 
This is the experimental signature that the taper becomes nearly single mode, as shown below. 

We have also observed with a video-microscope the exponential decrease of the taper waist, 
shown in solid squares in Fig |2] Note that the experimental points stops before the end of the 
tapering; indeed, a transverse motion of the thinned fiber in the turbulent microtorch gas flow, 
made the last image too fuzzy for a meaningful fiber diameter measurement. The experimental 
curve is well fitted by a exponential, as predicted by the model given in fSl, relying on volume 
conservation. According to this model, the taper shape in cylindrical coordinates, with origin 
z = at one end of the taper (see Fig. [TJ, is described by: 

r ro exp(-z//i) for 0<z<f 

r(z,L) = < , r , , (1) 

I w = ro exp {-L/2h) for I < z < 

with a symmetric profile for z > {L + h)/2. In this equation ro = 62.5 jJ-m is the initial fiber 
radius, L is the fiber lengthening and h the so called "hot zone" length. The "hot zone" is the 
cylindrical part where the softened fiber is drawn, roughly corresponding to the flame width. 
Though sensitive to the fiber location in the flame, the fitted h for adiabatic tapers remains in 
the range 6-8 mm, consistent with our flame geometry. 

Analysis of oscillations 

The observed oscillations originate from the beating of different local modes, which propagate 
along the taper with different propagation constants j3, =2n/X x A^eiT, i (^eff, i is their effective 
index), and recombine at the output with a relative phase which depends on the effective length 
of the taper. The amplitude of these oscillations varies during the lengthening of the fiber and 
reveals the efficiency of excitation and recombination of the different modes during the tapering 
process. In Fig|2] a peak-to-peak amplitude of 10% corresponds to an energy transfer of only 
5% from the fundamental to excited modes, resulting in a final transmission of 95%. 
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Fig. 3. Calculated effective index for the lowest modes of a cylindrical silica fiber, as a 
function of fiber radius. The thick solid and dashed lines are for HEj,, and HE2,, modes; 
the thin solid and dashed lines for the EHi„ and EH2,, modes, the thin dotted lines for the 
TEo„ modes. 



Figure [3] shows the effective indices of the first modes of a cylindrical waveguide with a 
uniform silica core and an air or vacuum cladding (A^ = 1). Due to strong guiding conditions, 
the relevant modes differ from the familiar LP modes, and include "hybrid modes" with more 
complex polarization properties. In particular, the fundamental mode is the hybrid HEn mode 
[|9l . For large radii the modes effective indices are close to the maximal value Ns = 1.453 
owing to a good confinement in silica and small numeric aperture. On the opposite for small 
radii the air penetration increases thus diminishing the effective index down to the minimal 
value min(A^eff) = 1, which is reached at mode cutoff. As a result, the effective index dispersion 
increases as the taper radius decreases, hence the gradual increase of the oscillation frequency. 

Further interpretation of the oscillations is obtained by analyzing the beating frequency of 
two given local modes, designated by indices / = 1 or 2. Note that due to the taper radius vari- 
ation along the propagation axis z, the associated propagation constants j3, (r) are z-dependent. 
Using the symmetry around the taper center, the accumulated relative phase writes: 

^niL)^2i^JJ Alin{r{z))dz + AjSnW^j ■ (2) 

A simple derivation taking into account the exponential variation of w vs. L, leads to the spatial 
angular frequency : 

^i2-^=A^i2(w)-^^(AM 



(3) 



This frequency can be derived by performing a short-time Fourier transform of the beating, 
known in acoustics as a sonogram or spectrogram 111 01 . Figure|4]shows such a spectrogram with 
the lengthening as abscissa and the normalized frequency KX/2n as ordinate. It should be em- 
phasized that a poorly adiabatic taper (with li <5 mm) has been intentionally chosen in order to 
make the discussed features more visible. In this figure is also plotted the normalized frequency 
calculated from eq. ^ and eq. ([T]i. They correspond to the beating of the fundamental mode 
HEii with the HE12, HE21 and TEqi modes, from the highest to the lowest curve, respectively. 
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Fig. 4. Transmittance curve and its short-time Fourier transform (window width is set to 
~ 0.5 mm). The solid curves are the frequencies calculated from eq. ^ with h as single 
fit parameter. The vertical lines underline the coincidence of amplitude drop with mode 
cutoff. 

Using the "hot zone" width h as an adjustable parameter, our model reproduces accurately the 
experimental observations, with the fitted value Iiq = 3.05 mm. 

As expected the most intense component arises from the beating of the fundamental mode 
with the HEi2 mode, the first excited mode that has the same symmetry (which is also the 
symmetry of the LPOl mode initially launched in the single mode fiber). The corresponding 
difference AA^eff is illustrated on Fig[3]by the righter-most arrow. Note that the small decrease 
of the spatial frequency just before cutoff arises from the inflection point of the dispersion 
curves visible in Fig. [3] and of the presence of the derivative term in eq. (O. 

After the cutoff of the HE 12 mode (occurring at w « 460 nm) the amplitude of the remaining 
oscillations is reduced to a few %, corresponding to the two lighter lines on the sonogram. They 
result from the weak excitation of two modes with a different symmetry, which can likely be 
attributed to imperfections in the early stage of tapering. When these two modes successively 
reach cutoff, the oscillations fully disappear. From then on the taper is singlemode, with a waist 
smaller than 300 nm. 

In order to confirm our interpretation of the oscillation pattern, we mapped the evanescent 
field of the taper in its final shape, similarly to the mapping experiment of the whispering 
gallery mode field reported in ifTTl . A thin fiber tip is scanned along the taper using a 3-axis 
piezo-translation stage while the intensity catched by the fiber tip is monitored. Looking only to 
the intensity distribution along the taper, we put the tip in contact with the taper This ensures a 
good mechanical stability and allows to scan over the whole 100 /im-stroke of the piezo stage. 

Figure m shows the detected signal as a function of position for three different places along 
the taper studied in Fig.|2l (1) in the central region, (2) at a distance 1.25 mm before the hot 
zone and (3) 2.25 mm before the hot zone. As expected, the average intensity becomes smaller 
when the taper becomes thicker, because the modes are more efficiently confined inside the 
silica. We also observe that the fringes visibility is much larger than on the taper transmittance 
curve: indeed the excited mode that carries only a few percent of the total power is less confined 
than the fundamental mode, and has therefore a larger relative intensity in its evanescent part. 
We performed a Fourier transform of these signals in order to identify the involved modes. The 
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Fig. 5. Near-field intensity mappings recorded (1) at tiie taper center , (2) 1.25 mm before 
the hot zone and (3) 2.25 mm before the hot zone . 

results are given in table [T] retaining in the detected signal and for each scan the two principal 
modes beating with the fundamental. The weights in table [T] give their relative contribution 
in intensity to the fiber tip signal. The periods A are recast in effective index differences by 
the relation ANf,ff = A/A, and the taper radius is determined by the co-existence of the two 
effective index differences taken into account (see black arrows in Fig.[3]l. 



Table 1. Frequency analysis of Fig.|5] 

Scan 1 Scan 2 Scan 3 

Offset position (mm) 1.25 mm 2.25 mm 

Weight 95% 4% 92% 3% 72% 20% 

AA^eff 0.22 0.16 0.35 0.30 0.26 0.40 

Beating mode HE21 TEOl HE 12 EH 11 HE 12 HE22 

Radius r(^m) 0.4 0.51 0.68 



Conclusion 

We have shown that the oscillations of the transmitted power during the tapering process pro- 
vides a real-time information on the modes propagating in the cylindrical part of the taper. The 
beating frequency is directly related to the associated effective indices. The successive drops 
observed on their intensity in the final stage is a convenient criterion to ascertain the cutoff of 
the last guided modes, and provides a mean to directly control the final taper diameter, in a size 
range where a direct optical measurement of the fiber diameter is no longer possible. Moreover 
this allows to control the fundamental mode effective index between 1.4 and 1.25, which is the 
appropriate range for microsphere and microtoroid's whispering gallery modes excitation. 
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